home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / soliton.ode.z / soliton.ode
Encoding:
Text File  |  1998-10-28  |  432 b   |  26 lines
  1. # You may run this example by doing:
  2. #
  3. #    ode < soliton.ode | graph -T X -C
  4. #
  5. # or alternatively, to create a real-time plot,
  6. #
  7. #    ode < soliton.ode | graph -T X -C -x 5 15 5 -y -3 3
  8.  
  9. # This example numerically solves an initial value problem arising in
  10. # John Selker's soliton thesis:
  11. # 
  12. # x'' = (2 / t) * x' + B * x
  13. # B = -1
  14. # x(5) = 0
  15. # x'(5) = 1
  16.  
  17. x' = vx
  18. vx' = 2 / t * vx + B * x
  19. B = -1
  20.  
  21. x = 0
  22. vx = 1
  23.  
  24. print t, x
  25. step 5, 15
  26.